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We point out that a pseudo-Newtonian interpretations of the (1/3) factor in the Sachs- Wolfe effect, 
which relates the fluctuations in temperature and potential, (ST/T) — (1/3)5<I>, is not supported by 
the General Relativistic analysis. Dividing the full gravitational effect into separate parts depends on 
the choice of time slicing (gauge) and there exist infinitely many different choices. More importantly, 
interpreting the parts as being due to the gravitational redshift and the time dilation is not justified 
in the rigorous relativistic perturbation theory. We suggest to regard the (1/3) factor as the General 
Relativistic result which applies in a restricted situation of adiabatic perturbation in the K = — A 
model with the last scattering occuring in the matter dominated era. For an isocurvature initial 
condition the corresponding result, (ST/T) — 2S^, has a different numerical coefficient. 

PACS number(s): 98.70.Vc, 98.80. Hw 



I. INTRODUCTION 

Some of the exact results in General Relativity (GR) 
can be given intuitive interpretations in terms of New- 
tonian concepts even though the structures of these two 
theories are very different, (i) The simplest example is 
the "derivation" of Schwarzschild radius by equating the 
Newtonian escape velocity to the velocity of light; one 
gets the correct numerical factor though the argument is 
not correct, (ii) The collapse time for a uniform density 
dust sphere under self gravity computed by Newtonian 
theory turns out to be exactly equal to the General Rela- 
tivistic result, (iii) The gravitational force on a material 
particle located inside the empty region of a spherical 
shell of matter vanishes in Newtonian theory. This is 
usually explained by the fact that the force in Newtonian 
gravity falls as (1/?"^) while the amount of matter inter- 
cepted by a a cone with fixed solid angle grows as (r^) 
thereby leading to cancellation of forces due to opposite 
pairs of material in the shell H . Incredibly enough this 
result is true in GR as well though the force law in GR 
is not strictly (1/r^) while the material intercepted by a 
fixed solid angle does increase as (r^). (iv) At a more 
subtle level, one can obtain the time-time component 
of Einstein's equations for FLRW (Friedmann-Lemaitre- 
Robertson- Walker) universe from Newtonian energy con- 
servation argument if we take the potential energy of a 
spherically symmetric mass distribution to be —{GM/R). 
The argument is again invalid in GR especially since 
Birkoff's theorem is applicable only for empty regions 
outside a mass distribution. In all the four cases men- 
tioned above, these are accidental coincidences and the 
heuristic arguments have no deeper significance. 

Recently, another pseudo-Newtonian argument has 
surfaced to explain the origin of (1/3) factor in the Sachs- 



Wolfe effect in cosmic microwave background (CMB) ra- 
diation H. The following heuristic argument is often 
seen in the literature including several textbooks [§|-^: 
(ST/T) = (l/3)<5$ = (5$ - (2/3)^$ where the first term 
arises from the gravitational redshift whereas the sec- 
ond term comes from the time dilation of the temper- 
ature perturbation ST/T = -5a/ a = -{2/3){St/t) ^ 
— (2/3)(S$ where a is the scale factor. For a similar argu- 
ment see [|0|. We point out in this paper that the above 
decomposition is motivated by a certain slicing condition 
(temporal gauge choice) out of infinitely many possible 
choices. More importantly, however, the interpretation 
of the terms as above cannot be explained even in that 
slicing or in any other slicing in the context of the cor- 
rect relativistic analyses, see below eq. (^). Hence, as the 
heuristic argument and interpretation are not supported 
by rigorous theory, we believe such are spurious and have 
at best the same status as the examples we mentioned in 
the first paragraph. 

Before proceeding with a rigorous analysis of the situ- 
ation, it is worth pointing out that there are some fun- 
damental difficulties in providing any Newtonian inter- 
pretation of the Sachs- Wolfe effect. Given any metric in 
GR, the Newtonian limit can be rigorously established if 
the metric can be reduced to the form goo = ~(1 + 2<I') 
and gap = Sap in some coordinate system. When one 
attempts to do this with FLRW metric, [see ||ll|l], one 
finds that the procedure is valid only for a region of size 
much smaller than the Hubble radius rf/f (t) = (d/a)^^. 
This merely reiterates the fact that at scales bigger than 
the Hubble radius one requires the full machinary of GR 
and - in particular - one needs to grapple with issues of 
gauge. Since the Sachs- Wolfe effect arises from scales 
which are bigger than Hubble radius at the epoch of re- 
combination one can anticipate that any Newtonian in- 



terpretation could have problem. This conclusion is also 
supported from the fact that one can consistently choose 
a gauge in perturbed FLRW universe in which 500 = ^ 1 • 
There is no Newtonian potential available in this gauge. 
The key reason for the above pseudo-derivation to work 
is because one introduced an ill-defined quantity {5a/ a) 
and linked it with 5T /T at one end [by Ta — constant] 
and with 5t/t at the other end [by a oc t^^^]. There is no 
rigorous interpretation of this quantity {5a/ a) possible in 
Newtonian gravity, or for that matter in the gauge invari- 
ant treatments of FLRW perturbations. In handling the 
cosmological perturbation one can adopt a heuristic ar- 
gument using 5a. However, in this way, one often ends up 
with an equation which is wrong (especially in a medium 
with pressure)|j. This is not surprising because perturb- 
ing the background system necessarily loses some degrees 
of freedom compared with perturbing the full system. 



II. TEMPERATURE ANISOTROPY IN AN 
ARBITRARY GAUGE 

We begin by obtaining the expression for the temper- 
ature anisotropy without imposing any choice of gauge 
so that we can study the results in the most general set- 
ting. We consider a spatially homogeneous and isotropic 
metric with the scalar-type perturbations 



ds^ = -a^{l + 2a)d7f - 2a^l3^adr]dx°' 



+a' 



gg](l + 2(^) + 27,„|^ 



dx^dxl^. 



(1) 



This represents a fairly general perturbed metric and no 
specific gauge has been chosen. The variables a, /?, if and 
7 are spacetime dependent scalar-type metric perturba- 
tions, and a vertical bar indicates a covariant derivative 
based on the background three-space comoving metric 

9ai3- ^^ introduce the variable x = o.{P + ^t) which 
gives the shear of the normal- hypersurf ace. (p is propor- 
tional to the perturbed curvature of the hypersurface, 
and a is the perturbed lapse function; we use notations 
in [0,|l3]. The combination x and the rest of the vari- 
ables used in the following are spatially gauge-invariant 
|l3| . The variables depend, however, on the temporal 
gauge (coordinate) transformation which corresponds to 
choosing the spatial hypersurface, i.e., the time slicing. 
Thus we have the freedom to impose a temporal gauge 
(slicing) condition which could be used as an advantage 
to handle the problems conveniently. The prime and the 
overdot indicate the time derivative based on 77 and i, 
respectively, with dt = adrj. 



The most general expression of the Sachs- Wolfe effect 
from the scalar-type perturbation can be found in eq. 
(15) of [Q. Ignoring the Doppler effects due to the ob- 
server's motion, and the emitting event (which is sub- 
dominant at large angular scales) the observable temper- 
ature anisotropy becomes 



5T 



5T^ 



(-X + -^) 



+ / {ax-fxYdy, 



(2) 



where the integration is along the photon's null geodesic 
path from the emitted {E) epoch to the observed (O) 
epoch. The variables a^ = a — x, 5T^ = 5T -t- HTy, 
and Lp^ ^ Lp — Hx are gauge-invariant combination^. 
tpy. becomes ip in the zero-shear gauge (often called the 
Newtonian gauge) which sets x = 0, etc. Using Bardeen's 
notation in [|l5| we have a^ = $yi and (p^ = ^h, [0- 
The gauge invariance of a combination assures that there 
remains no gauge (coordinate) mode and that its value 
remains the same in any gauge. However, it does not 
guarantee a certian gauge-invariant variable has an asso- 
ciated intrinsic physical meaning independent of the slic- 
ing condition. For example, there exist several variables 
ip^, ipsT = 'P + 5T/T, ip^ = ip — {aH/k)v where z; is a ve- 
locity related variable and fc is a wave number, etc., which 
are all gauge invariant. (An exception is Lpa based on the 
synchronous gauge which fixes a = 0). They all reduce 
to (p under the corresponding gauge condition which sets 
the variable in the subscript equal to zero. Although all 
these gauge invariant variables (we can make infinitely 
many different combinations) are the curvature variables 
in different time slicing (temporal gauge) we can safely 
regard them as completely different variables. 

Equation (g) was derived from the geodesic equations 
in the spacetime of eq. (0) without fixing the gauge condi- 
tion and without using the gravitational field equations. 
In the literature, the two terms in the RHS are often 
called the Sachs- Wolfe effect and the integrated Sachs- 
Wolfe effect, respectively. At this point it may be ap- 
propriate to quote a comment in [H which looks presient 
in the context of main point of this paper: "We empha- 
size again that in a generic gravitational field one cannot 
distinguish gravitational redshifts from Doppler shifts by 
any standard recipe; thus our division of equation . . . has 
only a heuristic significance^^ 

We shall now assume that: (i) the anisotropic stress 
can be ignored, so that we have[j a^ = —'Px'^ (n) I'^ = 
= A, and (iii) the medium is an ideal fluid with con- 
stant w = p//i (where p is the pressure and /i is the 
energy density) so that the growing solution of p^ re- 
mains constant in timq^. Given all these assumptions 



* This is in contrast with the perturbed Hubble parameter 
6H which has a rigorous geometric and kinematic meanings, 
and in fact is quite useful to handle the cosmological pertur- 
bation in a heuristic looking but fully rigorous manner |l2|l. 



^ For the gauge transformation properties, see eq. (2) in 

^See eq. (8) in ^ 
§ See eq. (18) in 



the integrated Sachs- Wolfe term vanishes, so that: 



6T 
~T 



^ + 2Hx 



4*^(7) 



(3) 



where we used ST/T\e — (1/4)5(^)|£; with 5(^) = 
(5/i(^)//i(^) denoting fractional energy-density perturba- 
tion of the photons. Now, the RHS is written without fix- 
ing the gauge yet, thus in a sort of gauge-ready form, but 
the sum is gauge-invariant. In this form we could under- 
stand why and how such a decomposition into the intrin- 
sic temperature perturbation, the gravitational redshift 
(or the time dilation) is dependent on the temporal slic- 
ing (gauge) condition of the spacetime. Although ip and 
X have meaning as perturbed curvature and shear of the 
normal three-hypersurface, and S(^j^ looks like a energy- 
density perturbation of photons, these variables acquire 
such a meaning only after fixing the temporal gauge (time 
slicing) condition where we have infinite choices. As men- 
tioned before the same variable evaluated in a different 
slicing (gauge) condition in general behaves as a com- 
pletely different variable. 

Equation (0) can be written in a suggestive form aa*j 



5T 



Vx 



aH 1^ 



(4) 



where v^ = v — {k/a)x and ^(-y)^ = S^^^ + A{aH/k)v. 
5tj\y is the same as Sij\ in the comoving gauge which 
sets the velocity variable v = 0; in a pressureless matter 
the test particles follow geodesies and hence the comoving 
gauge is equivalent to the synchronous gauge which fixes 
a = 0. Thus, in this form the variables are viewed (eval- 
uated) in mixed slicing (gauge). We can show that each 
of these gauge-invariant variables most closely resembles 
the Newtonian counterparts as: —fx^ '^i" ^^'^ ^x most 
closely reproduce the behavior of the perturbed gravita- 
tional potential i54>, the perturbed density Sn = Sp/p, 
and the perturbed velocity Sv in the Newtonian context, 
0,[|jl|. 

It should be obvious from the above two equations and 
discussion that the actual form of the terms in the right 
hand side depends very much on the gauge. It is best not 
to yield to the temptation of interpreting the individual 
terms "physically" - let alone try to fix the numerical 
prefactors. But if one insists on doing so, then the most 
natural choice is to interpret the first term --'/'x ^^ '^^^ 
to gravitational redshift (we have —ip^ = a^ — S^), the 
second term as due to Doppler effect and the third as 
arising from radiation field. As we shall see in the next 
subsection, even this interpretation is fraught with dan- 
ger but at least the coefficient of —fx is now unity. 



III. ADIABATIC PERTURBATIONS 

We consider a system with radiation (7) and matter 
(to). The adiabatic condition 



S = (5(„i) - -(5(^) = 0, 



implies 



c _ Sp 1+ R 3p{ni) 

^ — — ": , „ '« "(7)i ^ ■ 



p 1 + AR/i 



4^(7) 



(5) 



(6) 



thus we havq^ 5 



{j)v 



,j = ^{k/aH)'^(px which is sub- 
dominant in large angular scales corresponding to the 
large-scale k/aH <^1 at E. We have|^ 



2 1 k 



(7) 



for the growing solution. Thus, adding the first two terms 
in the RHS of eq. (0) we finally have 



^y(Ad) 



T 



l + 3w 



3(1 



r'/'x 



-3^x 



. - 3^^' 



(8) 



where the second equality follows by assuming matter 
domination at E, so that w = Opi. In our case S^ does 
not evolve in time. We stress again that the result in 
eq. (0) is valid under many conditions mentioned above, 
especially the ones above eq. (^); e.g., the simple result 
in eq. (pi) does not hold (for example) in a model with 
additional cosmological constant, and in such a case we 
should go back to the general form in eq. (||). 

It was often stressed in the literature that the —(1/3) 
factor comes directly from the metric part in the syn- 
chronous gauge whereas it gets a contribution of — 1 from 
the metric and (2/3) from the intrinsic temperature part 
in the zero-shear gauge [|l^,D . The origin of such gauge- 
dependent interpretations can be understood simply by 
rewriting the RHS of eq. (pi) or eq. (0) in the respective 
gauge conditions. In the zero-shear gauge, x = 0, we 
have 



ST 
~T 



-"Px 



h. 



(7)X- 



(9) 



Notice that in the large-scale the temperature part 5i^\^ 
is dominated [when viewed in the comoving gauge, com- 
pare with eq. (0)] by the metric, and does not behave 
like an ordinary temperature. Instead, it gives |</3xi and 



** With the wave number k appearing in the equation the 
variables can be regarded as the Fourier transformed ones. 
To the linear order the same equations in configuration space 
remain valid in Fourier space as well. Thus, we ignore specific 
symbols distinguishing the variables in the two spaces. 



ttSee eq. (6) in & 
** See eq. (7) in^. 

^^ Using ^ = 3(i+m) ^ t'^^ heuristic argument mentioned 
in the introduction also produces the result for general w at 



the emission epoch E H. 



we use —1 + 1 = — i to get the final result. This is a 
rigorous argument, and one should not confuse this with 
the heuristic one mentioned in the introduction; except 
for the similar division into —1 and | the origins and the 
interpretations are completely different. Therefore, the 
heuristic interpretation is not based on this zero-shear 
gauge analysis. The synchronous gauge coincides with 
the comoving gauge in the matter dominated era (MDE), 
thus in the comoving gauge we have 



ST 



1 



(-(/j^, + 2Hxv) + -5. 



(l)v- 



(10) 



ideal fluid with w — constant we have ip^ — 
for the growing mode. The temperature part 



For an 

3+3u) "^X 

now behaves like the conventionally known temperature 
fluctuation, thus is negligible in the large angular scales 
compared with the potential fluctuation. Therefore, us- 
ing eq. (0) with v^ — —{k/a)xv the metric part gives 
— |(Py directly. For the original derivation see 0; see 
also 



IV. ISOCURVATURE CASE 



The isocurvature condition is 5^v = under which 
we have Lp-^ — Q and Lp.^ = 0. This condition implies 
S = — |(i?~^-(-l)(5(^)i,. The isocurvature initial condition 
is imposed early in the radiation dominated era (RDE). 
Einstein's equations givef"] ip^ — —H{p.+p)^^Spv which 
shows that the initial isocurvature perturbation can gen- 
erate ify. For an isocurvature mode we have Spy = 
— (l/3)^(m)(l + R)^^S. Assuming the last scattering 
epoch E occurred in the MDE we have 



_ 1 /-^DE dR _ 1 



(11) 



at E where we used S — constant in the large-scale limit 
pl[ ; this argument was used by Liddle and Lyth in [g2| . 
Thus, in the MDE we have ip-^ = ^pv = ^S; this shows 
the amount of curvature perturbation (p-^ in MDE gener- 
ated from the initial isocurvature perturbation S in RDE. 
In MDE we have jS(^^^y = —^S. Therefore, from eq. (||), 
using eq. (0) and assuming MDE at E, we have 



^y(Iso) 



T 



'"^V, 



2<5$, 



fl2) 



which is six times larger than the adiabatic result. For 
original derivations, see below eq. (3.5) of [g3| and below 
eq. (5.27) of pa. 



V. DISCUSSIONS 

It is clear from cq. (ph , that we can divide the terms in 
different ways (which is actually what the gauge choice is 
doing). Correspondingly there are many different ways to 
reach t/ie same final results in eqs. (@,0). In other words, 
we have infinitely many different ways of introducing slic- 
ing, thus viewing each variable in different gauges. While 
doing an actual calculation we need to choose the gauge 
(as we mentioned, a gauge- invariant variable is equivalent 
to a variable based on a certain slicing condition which 
fixes the gauge mode completely), but the final physical 
results should be the same independently of which gauge 
we have chosen. Our results in eqs. (p|,p^) are the fi- 
nal results where 5$ can be interpreted as the perturbed 
Newtonian potential which is related to the density con- 
trast through Poisson's equation. 

The pseudo-Newtonian method described in the intro- 
duction is closely related to the decomposition in cq. (||). 
However, as we have shown below eq. y^ such an in- 
terpretation is not supported by analyses in that gauge, 
which is true even in the context of our gauge-ready form 
in eq. (||). That is, interpreting the parts as being due 
to the gravitational redshift and the time dilation is not 
justified in the rigorous rclativistic perturbation theory. 
Hence it is difficult to imagine that such a heuristic argu- 
ment captures the basic physics. We believe it is yet an- 
other curious coincidence in General Relativity in which 
a pseudo-Newtonian argument does lead to the correct 
final result. 

In a classic book by Zel'dovich and Novikov [^ we 
find a statement: "ffowever, . . . the gravitational shift 
contains the factor 1/3; it is still unclear how to interpret 
this coefEcient classically." which still seems to be true. 
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*** This follows from eq. (18) in |l4| and the conservation 
property of pv under an adiabatic condition: p^ ~ C. Or, 
see eqs. (50,51) in ||. 

tt+Tliis follows from eqs. (12,14,18) in lli. 
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